Testing the itinerancy of spin dynamics in superconducting Bi2Sr 2 CaCu 2 08+(5 
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Much of what we know about the electronic 
states of high-tempe ratur e superconductors is 
due to photoemissiorP^II] anc j scanning tunnel- 
ing spectroscopy 4 5 studies of the compound 
Bi2Sr2CaCu20g +d -. The demonstration of well- 
defined quasiparticles in the superconducting 
state has encouraged many theorists to apply 
the conventional theory of metals, Fermi-liquid 
theory, to the cuprates. 6 * 7 * 8 * 9 ! In particular, the 
spin excitations observed by neutron scatter- 
ing at energies below twice the superconduct- 
ing gap energy are commonly believed to cor- 
respond to an excito nic state involving itinerant 
electrons J 10 * 11 * 12 * 13 * 14 ! Here, we present the first 
measurements of the magnetic spectral weight 
of optimally-doped Bi 2 Sr2CaCu 2 8+ 5 in absolute 
units. The lack of temperature dependence of 
the local spin susceptibility across the super- 
conducting transition temperature, T c , is incom- 
patible with the itinerant calculations. Alter- 
natively, the magnetic excitations could be due 
to local moments, as the magnetic spectrum is 
similar to that in Lai .875Bao.i 2 5Cu04p^ where 
quasiparticles 16 and local moment^ coexist. 

Bi 2 Sr2CaCu208+(5 has been the cuprate system of 
choice for surface-sensitive techniques such as angle- 
resolved photoemission spectroscopy (ARPES) and scan- 
ning tunneling spectroscopy (STS) because it cleaves eas- 
ily, thus allowing simple preparation of fresh surfaces. 
ARPE^MUand STp studies of Bi2Sr 2 CaCu 2 8+l 5 have 
convincingly demonstrated the existence of coherent elec- 
tronic excitations (Bogoliubov quasiparticles) in the su- 
perconducting state. The excitation gap for quasiparti- 
cles has d-wave symmetry, going to zero at four nodal 
points along the nominal Fermi surface in the t wo- 
dimensional reciprocal space for a Cu0 2 plane On 
warming into the normal state, coherent electronic states 
are observed, at most, only over finite arcs about the 
nodal points; in the "antinodal" regions, there is a, so- 
called pseudogap and an absence of quasiparticles 

The ARPES results on Bi 2 Sr2CaCu208+5 have been 
used as a basis for predicting the magnetic excitation 
spectrum in the supercond uctin g state, measurable by 
inelastic neutron scattering! 10 ! 14 ! To detect the magnetic 
signal with neutrons, however, one needs crystals of large 
volume in order to compensate for limited neutron source 
strength and weak scattering cross section, and such crys- 



tals have been difficult to grow. The magnetic spectral 
weight is typically presented as the imaginary part of 
the dynamical spin susceptibility, x"(Q,w); here Q is 
the wave vector of a magnetic excitation and E = two is 
its energy, with fi being Planck's constant divided by 2tt 
and u> the angular fr equency. Previous neutron scatter- 
ing studie d 11 ' 12 1 13 - 11 ^ of Bi 2 Sr 2 CaCu 2 8+5 , working with 
crystals of limited size, focused on the change in x"(Q; w ) 
on cooling through the superconducting transition tem- 
perature, T c . 

A sustained effort has finally yielded crystals of suffi- 
cient size^to allow a direct measurement of x"(Q, w) on 
an absolute scale. The as-grown crystals correspond to 
the condition of "optimal" doping (maximum T c ), with 
T c = 91 K. To des cribe our results, we will follow previous 
practice^ 11 ! 12 ! 13 ! 14 ! and make use of a pseudo-tetragonal 
unit cell, with lattice parameters a = 3.82 A (parallel to 
in-plane Cu-0 bonds) and c = 30.8 A. Wave vectors Q 
are given in reciprocal lattice units of (2n/a, 2n/a, 2ir/c). 
The antifcrromagnetic wave vector within a Cu0 2 plane 
corresponds to Qaf = We note that the crys- 

tallographic unit cell is larger and orthorhombic, with a 
long-period incommensurate modulation in one in-plane 
direction. The substantial atomic displacements asso- 
ciated with the incommensurability can modulate elec- 
tronic and magnetic interactions^ however, we did not 
detect any anisotropy of the magnetic scattering associ- 
ated with the incommensurate modulation direction. 

The magnetic response was measured by inelastic neu- 
tron scattering using a time-of-flight (TOF) technique 
(see Methods). Figure [T] presents constant-energy slices 
of the scattered signal as a function of momentum trans- 
fer about Qaf for several energies and temperatures. 
The results are plotted in the form of x"(Q,^), as dis- 
cussed in Methods. Besides the magnetic scattering, 
peaked at Qaf, the measurements also include contribu- 
tions from phonons; the phonon contribution should be 
approximately independent of temperature, but should 
vary, on average, as Q 2 . In this figure, we have subtracted 
a fitted background of the form C0 + C2Q 2 , where the coef- 
ficients Co and C2 depend on energy, with C2 independent 
of temperature. (The magnitude of the background can 
be seen in Fig. [2^-c, where plots of the measured inten- 
sity along specific directions in Q are plotted, with the 
fitted background indicated by a dashed line.) The mag- 
netic signal is expected to be weaker and more diffuse at 
310 K, and the signal found there is generally consistent 
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FIG. 1: Constant energy slices through the scattered 
signal. Time-of-flight data (for Ei = 120 meV) converted to 
the form of x"(Q> w ) ( see Methods), with background sub- 
tracted (see text). Data from four equivalent zones have been 
averaged to improve statistics. Each column corresponds to 
a different temperature: 10 K, 100 K, and 310 K, from left to 
right. Each row shows a different energy transfer: 36 meV, 
42 meV, 54 meV, and 66 meV, from bottom to top. For each 
slice, the data have been averaged over an energy range of ±3 
meV. White indicates areas not covered by detectors. 



with this expectation. The diagonal streak of scattering 
for hu — 36 meV at 310 K indicates a phonon contribu- 
tion with intensity modulated in Q. 

For a more quantitative analysis of the magnetic sig- 
nal and its temperature dependence, we can take cuts 
through these images and fit gaussian peaks to the struc- 
tures. Examples of such plots are shown in Fig. [2^-c 
and Fig. [3] From the fits, we obtain both the magni- 
tude and Q dependence of the magnetic response. The 
effective dispersion of the magnetic excitations is plot- 
ted in Fig. [2J3. Within the uncertainties, it is effec- 
tively isotropic about Qaf- At low temperature, the 
overall shape of the dispersion is qualitatively consistent 
with the hour-glass spectrum that appears to be univer- 
sal among the cupratesP^ with the excitations crossing 
Qaf at hu> ~ 40 meV^ the response falls off below a 
spin gap energy of ~ 25 meV (see Fig. pjl. The effec- 
tive upward dispersion, above 40 meV, is comparable to 
that observed^ in underdoped YBa 2 Cu306.5, and much 
steeper than thaP^in optimally-doped YBa2Cu3 06.95- 

The method that we have used to identify the magnetic 
response and separate it from the large phonon back- 
ground is dependent on assumptions. To get a direct 
measure of the magnetic response, we have performed 
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FIG. 2: Analysis of magnetic response, a-c, Cuts 
through x"(<i> (for Ei — 120 meV) along q = Q — Qaf, av- 
eraged over q along [1, 0, 0] and [0, 1, 0] (paths A and C, shown 
in d). Energy transfers are: a, 66 meV; b, 54 meV; c, 42 meV. 
Black symbols: T = 10 K; red symbols: T = 100 K; error bars 
indicate standard deviations. The solid lines through the data 
represent fitted gaussian peaks — either a single peak or a pair 
of peaks constrained to be symmetric about q = 0. Dashed 
lines indicate background, e, Plot of peak positions obtained 
from the peak fitting. Circles: Ei — 120 meV; diamonds: 
Ei — 200 meV; bars indicate peak widths. For Tlui > 40 meV, 
points represent fits to the averaged A and C cuts; virtu- 
ally identical results were obtained for B and D cuts. For 
hu> < 40 meV, only fits to D cuts (as in Fig. [3f were used, be- 
cause of the more complicated phonon background. Solid blue 
lines indicate magnetic dispersion for YBa2Cu3 06.95 from 
Reznik et aZp^ gray lines indicate results for YBa2Cus06.5 
from Stock et aZ™ 



a second experiment with spin polarized neutrons (see 
Methods and Supplementary Information). Because of 
intensity limitations associated with this technique, our 
useful results arc limited to measurements at Qaf ; these 
are presented in Fig. |4ji. The lower panels compare re- 
sults for x"(Qaf, extracted from fits to the TOF data. 
The data in Fig. [4^, and c are generally consistent, with 
a substantial temperature-dependent change in signal at 
fvjj ~ 40 meV. More modest changes with temperature 
are seen in Fig. [4j), which corresponds to measurements 
with a different incident neutron energy, E{. 

To understand the sensitivity to Ei, we 
have to consider the fact that Cu02 planes in 
Bi 2 Sr 2 CaCu208+5 come in correlated bilayers, just 
as in YBa2Cu306+a;. The magnetic response from 
a pair of layers can be separated into components 
with symmetry that is odd or even with respect to 
exchange of the layers. The relative cross sections for 
the odd and even responses depend on the momentum 
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FIG. 3: Magnetic response at lower energies Cuts 
through x"(q,w) along q = [1,-1,0] (path D in Fig. |5Ji) 
for excitation energies of 24, 30, 36, and 42 meV (bottom to 
top) measured at T — 10 K with Ei = 120 meV. Data sets 
have been offset vertically for clarity. Dashed lines indicate 
background. Error bars indicate standard deviations. 



transfer component along the c-axis, Q c . For the TOF 
measurements, Q c varies with hto in a fashion that is 
determined by the choice of The relative weight 

of the odd contribution as a function of £7, and tuv is 
indicated by the dashed gray lines in Fig. [4] the sum 
of the weights for the odd and even components at any 
given Tiuj is equal to one. For JSj = 200 meV (Fig. Hp), 
the response at hui ~ 40 meV is mostly odd. For the 
polarized-neutron study, all of the measurements were 
done at a fixed wave vector that maximizes the odd 
component. For E± — 120 meV (Fig. ^p), there is a 
bigger contribution from the even component, which 
shows less change with temperature.^ 

The enhancement of the odd component of x"(Qaf, w) 
in the superconducting state for hui ~ 40 meV is consis- 
tent with previous worl^J^l however, our absolute mea- 
surements of the magnetic response also indicate signifi- 
cant weight in the normal state. To get another perspec- 
tive on the temperature dependence, we have integrated 
the fits to the magnetic response over the in-plane mo- 
mentum transfers (assuming %"(Q,oj) to be isotropic in 
q = Q Qaf) to yield the local susceptibility, x"{uo), 
plotted in Fig. HJi-f. As one can see, there are minimal 
changes in x" (ui) between 10 K and 100 K. Thus, the 
changes in x"(Qaf, w ) appear to be associated with a 
narrowing of the Q dependence at low temperature. 

Our results cast considerable doubt on the common 
view of the ~ 40-meV "resonance" mode as a spin-1 ex- 
citation of a pair of quasiparticles in the superconduct- 
ing statep^ where the excitation is presumed to be res- 
onant or excitonic because it occurs at an energy below 
twice the maximum superconducting c?-wave gap energy 
(2A). From ARPEp and STp studies, we know that 




20 40 60 80 20 40 60 80 100 



h w (meV) 

FIG. 4: Temperature and energy dependence of x" • 

a— c, ^"(Qi 1 ^) a l Q = Qaf, with obtained from: a, spin- 
polarized beam measurements at Q = (0.5,1.5,4.5); b, TOF 
measurements with Ei = 120 meV; c, TOF data with Ei = 
200 meV. Error bars indicate standard deviations, d— f, Lo- 
cal dynamic susceptibilty, x"(uo), for: d, Ei = 80 meV; e, 
Ei — 120 meV; f, Ei — 200 meV. For all panels, black circles: 
T = 10 K; red diamonds: T = 100 K; vertical bars indicate 
standard deviations, based on the least-squares fits used to 
evaluate and integrate the Q dependence of Q), while 

horizontal bars indicate energy range over which the signal 
was averaged. Inset in d shows a rough estimate of 
at 10 K for odd (dashed) and even (dotted) components; at 
100 K, gap is decreased, widths increased, areas kept approx- 
imately constant. These estimates are used (with amplitude 
adjustments for a— c) to produce the solid lines in each panel, 
taking into account the weight variations of the odd and even 
components with hui and Ei\ gray dashed lines show the rel- 
ative weight of the odd component (magnitude oscillates be- 
tween and 1). 



2A w 80 meV for optimally-doped Bi2Sr2CaCu20s+5. It 
has been established that the superconducting gap de- 
creases to zero at T — T c for wave vectors on the nodal 
arcP A pseudogap remains at antinodal wave vectors 
for T > T c , but this involves only a weak depression 
of the electronic spectral function at the Fermi energyP 
Thus, if the magnetic excitations were due to quasiparti- 
cles, we should have seen dramatic changes in x"( w ) be- 
tween 10 K and 100 K over the entire measured energy 
ranged clearly, our results are inconsistent with such 
a picture. Similar inconsistencies for an electron-doped 
cuprate have been emphasized by Kriiger et oZPS Fur- 
thermore, the strong temperature dependence of x"( w ) 
inferred from a conventional quasiparticle-based analy- 
sis of optical conductivity dataP is inconsistent with our 
direct measurement, raising a challenge to that analysis. 

The magnitude of x"^) m the energy range 
of 40-50 meV (Fig. g}i-f) is smaller than, but 
within a factor of two of, that for YBa2Cu.3O6.95pS 
Lai.84Sr . 16 CuO4pS! and LaL875Bao.125CuO4.ES In the 
case of La2- a; Ba a ;Cu04, it has been shown that the mag- 
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netism is dominated by local moments on Cu atomsf^ 
Given the comparable magnitude of the magnetic re- 
sponse, universal dispersion, and the limited sensitivity 
to the presence of a superconducting gap, it seems likely 
that the spin fluctuations in I^Si^CaC^Os+a are due 
primarily to local-moment magnetism. A recent study of 
YBa 2 Cu3 06.95 reached a similar conclusion.^ 

Lai.875Bao.i25Cu04 may seem an unlikely case to com- 
pare with optimally-doped Bi 2 Sr2CaCu208+,5, as th e for- 
mer compound exhibits charge and spin stripe order p2H2] 
Nevertheless, a recent ARPES studjBSI has demonstrated 
the coexistence of nodal quasiparticles with the localized 
electronic moments of the spin stripesp^ and there is ev- 
idence for two-dimensional superconducting correlations, 
as well. 28 Explaining such counterintuitive behavior re- 
mains a challenge for theorists. 

Methods 

The Bi 2 Sr2CaCu208+(5 crystals were grown at 
Brookhaven using the traveling-solvent floating-zone 
methodP Plate-like single crystals with sizes up to 
50 x 7.2 x 7 mm 3 were obtained. Magnetic susceptibil- 
ity measurements indicate an onset of diamagnetism at 
T c — 91 K. Four crystals, with a total mass of 19 g, were 
mounted on aluminum supports and co- aligned for the 
experiment on the MAPS time-of-flight spectrometer at 
the ISIS spallation facility, Rutherford Appleton Labora- 
tory. The sample was mounted in a closed-cycle He refrig- 
erator for temperature control, with the c-axis oriented 
along the incident beam direction. The combinations of 



(incident neutron energy Ei, Fermi chopper frequency, 
typical integrated beam current) used for the measure- 
ments were (80 meV, 250 Hz, 5000 fiA h), (120 meV, 
350 Hz, 9000 fiA h), and (200 meV, 500 Hz, 7000 ^A 
h). Measured differential cross sections, d 2 a / dH.dE , were 
converted to absolute units by normalization to measure- 
ments on standard vanadium foils. The dynamic suscep- 
tibility was extracted using the formula 

where /cu(Q) is the anisotropic magnetic form factor for 
CuP 

A second set of 7 crystals, with a mass of 18.75 g, 
was aligned with x-rays at Brookhaven (collective mo- 
saic < 2°) for the experiment on IN22 at the ILL. 
The spectrometer is equipped with a vertically-focusing 
monochromator and horizontally- focusing analyzer, both 
made of Heusler alloy crystals. The sample, oriented 
with [100] and [039] directions in the horizontal scatter- 
ing plane, was mounted in an ILL He- flow cryostat, inside 
of Cryopadum, a system providing spherical polarime- 
try capabilities.^ The neutron polarization flipping ra- 
tio measured on a strong Bragg peak was 16. To obtain 
the inelastic magnetic scattering cross section, we used an 
appropriate combination of spin-flip intensities measured 
for three orthogonal neutron polarizations (see Supple- 
mentary Information) with a fixed final energy of 30.5 
meV. The results were normalized to the time-of-flight 
data through non-spin-flip measurements of phonons. 
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SUPPLEMENTARY INFORMATION 
Polarization analysis 

Let Q (in the horizontal scattering plane) define the x 
direction, and let z be vertical; then y is the orthogonal 
direction in the scattering plane. Using CryopadpEl it is 
possible to set the initial or final polarization direction to 
be along any of these directions. (Note that the sample 
sits in a zero-field, magnetically-screened environment.) 
For our problem, it is sufficient to consider only cases 
where the initial and final polarizations are parallel. 

We can measure scattered intensities as a function of 
polarization orientation (xx, yy, or zz). If the flipper 
is off, then we will detect neutrons that have the same 
spin orientation as in the incident beam (designated ++ 

or ); with the flipper on, we detect spins that have 

been flipped by scattering from the sample (-1 — or — V). 
Thus, we can label the measured intensities as 1^%, where 
a — x, y, z, and s = +, — . 

For the scattering cross section from the sample, we 
will use (Jaa as an abbreviation for the differential cross 
section; note that it is labeled with the same polarization 
orientation and spin-direction indices. Let N represent 
the contribution due to nuclear scattering (including con- 
tributions from the sample holder, cryostat, etc.), and 
let M a represent the magnetic contribution due to fluc- 
tuations along the a direction. (Note that M x is zero 
because those fluctuations are parallel to Q, and only 
the fluctuations perpendicular to Q have a finite cross 
section.) For I^S^CaC^Os+a, there should be no sig- 
nificant nuclear spin incoherent scattering, so we will as- 
sume that the background B depends on counting time, 
but not on spin direction. With these definitions, all 
of the relevant cross sections can be written in a simple 
from: 
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From the spin-flip (SF) cross sections alone, we can 
uniquely extract the magnetic cross section with the com- 
bination 

M y + M z = 2a+- - a+- - a+7 ■ (8) 

Of course, the measured intensities I^ s a are limited by 
the degree of polarization achieved and maintained by 
the monochromator, analyzer, and the guide fields. (We 
will assume that the spin flipper performs perfectly.) To 
collectively characterize the effective polarizing efficiency, 
one typically measures the flipping ratio R by sitting on a 



Bragg peak and measuring the ratio of the intensity with 
the flipper off to that with the flipper on (say, Z+, + / I+T ) . 
By evaluating the measured intensities, taking account of 
the finite beam polarizations, one finds that 

My + M z = |±| (2I+- - I+- - I+-) . (9) 

This is the formula that we used to extract the magnetic 
cross section from measured SF intensities for the three 
polarization directions. 



